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ABSTRACT: In determining the natural modes and jreguencies oj a lin- 
ear elastic structure, Gut/an reduction is often used to reduce tne size of 
the mass and stljjness matrices and tne solution of the reduced system is 
ootalned Jlrst. The reduced system modes are then expanded to toe size oj 
tne original system by using a static transjormation itniclng tne retained 
degrees oj jreedom to tne omitted degrees of freedom, in the present paper, 
tne transformation matrix oj Guyan reduction is modljied to include ad- 
ditional terms from a senes accounting jot tne inertial ejjects. However, 
the inertial terms are dependent on tne unlcnoum frequencies. A practi- 
cal approximation is employed to compute the inertial terms urttnout any 
iteration. This new transjormation is implemented in NASTRAN using 
a DMAP sequence alter. Numerical examples using a cantilever beam il- 
lustrate tne necessary condition Jor allotting a large number oj additional 
terms in the proposed senes correction oj Guyan reduction. A practical 
example oj a large model oj the Plasma Motor Generator module to be 
jiou m on a Delta launen vehicle is also presented. 

1. Introduction: The dynamic analysis of complicated structures often produces large 
finite element models. In some instances, the automated computer procedures to generate 
finite element meshes also lead to large models. These highly refined models are really a 
byproduct of the use of model generating software and they may not be needed for accuracy. 
A common approach to reduce the size of the eigenvalue problem for structural dynamics 
applications is Guyan reduction. This approximate method finds its place among other 
applications also. For the purposes of cost-effectiveness, Guyan reduction is employed 
in Coupled Loads Analysis using substructuring techniques. In the experimental modes 
analysis, analytical selection of retained degrees of freedom for Guyan reduction is used as 
a guide to select accelerometer locations on the test article. Mass weighted orthogonality 
computations between the test and analytical modeshapes are performed using Guyan 
reduction. 

While Guyan reduction [1] is exact in static applications, it introduces approximations 
in structural dynamics. The correct relationship between the retained and omitted degrees 
of freedom can be expressed in the form of a series. The Guyan reduced mass and stiffness 
matrices, available in explicit form, are used to compute the Beries terms approximately. 
The Guyan reduced matrices provide the best possible solution without requiring any 
further iterations. The condition for convergence of the series and the relationship of this 
series transformation to the improved reduced system (IRS) introduced by O’Callahan [2] 
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are examined in this paper. 

2. Theory: The eigenvalue problem from the structural dynamic analysis is given as 

Ku = XMu (1) 

Eq. (1) can be written in the partitioned form as, 
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where u„ represents the eigenvector of the retained degrees of freedom and v 0 the eigen- 
vector of the degrees of freedom omitted in the Guyan reduction. Mij and Kij are the 
corresponding submatrices of the mass and stiffness matrices respectively and A is the 
eigenvalue. The second partition of Eq. (2) can be written separately as 


(Koa — XM oa )u a + ( K 00 — XM 00 )u 0 — 0 
Expanding the vector u 0 in terms of u a from Eq. (3), 

U 0 = - ( Koo - XMooY 1 (Koa - XM oa )u a 

= - (/ - XK-'Moo )' 1 (K-'Koa - W-'Moa) 


(3) 
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Guyan reduction transformation leaves out the frequency dependent terms in Eq. (3). 
Hence, the regular Guyan reduction transformation becomes, 

= ~K 0 g K oa u a (5) 

If the condition for convergence (Section 4.) is satisfied, the inverse of (7 - \K~'M \ 
can be expanded in Neumann series as, °° °° 


(/ XK 00 l Moo) = I + XK~oM 00 + A 2 [K 0 * M ob ] 2 + . . . 

Using Eq. (6) in Eq. (4) and simplifying the terms yields, 

*o = - [K-'Koa + BX + ABX 2 + A*B\ 9 + ...]u a 

where 

A = K; 0 : M 00 and B = K~ l M oa - AK~ l K, 


l oo dX oa 


( 6 ) 
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The exact relationship between *. and u 0 in Eq. (7) involves nonlinear terms of the un- 
known eigenvalues (A). A practical approximation to compute these terms in Eq. (7) can 
be made from regular Guyan reduction by taking, V 


K r U a =S A M r U t 


(9) 
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where K r and M r are Guyan reduced stiffness and mass matrices respectively and are 
given explicitly as, 


Kr = K a a - KtoR-gKoa 

M r — M a a ~ M a oK 0 o Koa ~~ Koodoo * d" K*oK 0 o M 0 o^ 0 o ^oa 


( 10 ) 


From Eq. (9), it is seen that, 

A« a = M^KrU a (H) 

Using Eq. (11) repeatedly, it can be shown that, 

X’«. = (M-'Krf*. 

: ( 12 ) 

A’v a = C'u a , C = M~ l K r 

Substituting Eq. (12) into Eq. (7), the relationship between * in Eq. (1) and * a becomes, 

u = Tu a (13) 


where 


/ 


T = 


-K-'Ko o 


+ ^ A*~ l BC* 

«'= 1 . 2 ,.. 


(14) 


By applying the relation between u and * 0 in Eq. (13), the new improved matrices from 
series reduction can be obtained as, 


TZ = T t KT and H = T t MT 


(15) 


It is interesting to note that M oa vanishes for lumped formulations of the mass matrix. 
Taking the value of t to be unity, the transformation in Eq. (14) reduces to 


T = 


-K~}K oa -f BC 


which is the improved reduced system (IRS) proposed by O’Callahan [2]. 


( 16 ) 


3. DMAP Alter: A rigid format alter for dynamic analysis in NASTRAN has been 
developed to incorporate the improved Guyan reduction with the series terms. A parameter 
called GOPT is used to choose the number of correction terms. The alter listing is also 
provided in this section. 
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$ CSA/NASTRAN ALTER FOR IMPROVED GUYAN REDUCTION 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

$ 

RFINSERT SMP2 $ 

PARAM //C,N,NOP/V,Y,GOPT=-l $ 

PARAM //C,N,SUB/V,N,GOUT/V,Y,GOPT/C,N,2 $ 

COND LGOPT.GOPT $ 

UPARTN USET,MFF/MAAB,M0A„M00/*F7*A7*0* $ 

FBS LOO„MOO/AMAT/l $ 

FBS LOO„MOA/BMATl/l $ 

MPYAD AMAT,GO,BMATl/BMAT $ 

SOLVE MAA.KAA/CMAT/l $ 

$ 

MPYAD BMAT,CMAT,/SUM $ 

COND OUT, GOUT $ 

MATMOD SUM,„„/PRDT,/13 $ 

LABEL LOOPTOP $ 

EQUIV SUM, SUMl /NEVER $ 

EQUIV PRDT,PRDTX/NEVER S 
SMYPAD AMAT,PRDT,CMAT,„/PRDTX/3 $ 

ADD SUM,PRDTX/SUM1 $ 

EQUIV SUMl, SUM/ ALWAYS $ 

EQUIV PRDTX,PRDT/ALWAYS $ 

REPT LOOPTOP, GOUT $ 

LABEL OUT $ 

ADD GO, SUM/GONE $ 

SMP2 USET,GONE,MFF/MAA $ 

SMP2 USET,GONE,KFF/KAA $ 

LABEL LGOPT $ 


4. Validity of Guyan Reduction: The inverse of the matrix [I - X K~JM 00 ] in Eq. (4) 
can be expanded as a converging Neumann series only if all the eigenvalues of A K~„M O0 are 
less than unity. In other words, the Guyan reduction is valid only for those frequencies les6 
than the smallest frequency of the eigenvalue problem formed out of the omitted degrees 
of freedom. The effect of violating this condition will be scrutinized in the next section. 

5. Demonstration Examples: 

5.1 Uniform Cantilever: The first example is concerned with a cantilevered bar clamped 
at one end. The relevent structural parameters are taken to be the modulus of elasticity 
( E ) being equal to 30 x 10 8 psi, weight density ( pg ) as 0.2839 Ib/in 3 , area of cross section 
as 1 in 2 and the length of the bar ( L ) as 72 in. 

The characteristic equation of this cantilever i6 cos y/p/EuL = 0 from which the 
theoretical natural frequencies can be computed. The cantilever is divided into twenty finite 
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elements. The retained degrees of freedom for Guyan reduction are the axial displacements 
at the free end and at two successive nodes. The reduction transformation includes n as 
the number of additional series terms. The natural frequencies from improved Guyan 
reduction for different values of n are listed in Tables 1 through 5. 


Tfeble 1. Cantilever Frequency Comparisons (n = 0) 


Standard Guyan Reduction 


Mode 

Number 

Theoretical 
Frequency (Hz) 

Computed 
frequency (Hz) 

Error 

% 

1 

7.012E2 

7.428E2 

5.926E0 

2 

2.104E3 

7.562E3 

2.595E2 

3 

3.506E3 

1.655E4 

3.722E2 


Table 2. Cantilever Frequency Comparisons (n = 1) 


Mode 

Theoretical 

Computed 

Error 

Number 

Frequency (Hz) 

Frequency (Hz) 

% 

1 

7.012E2 

7.012E2 

8.014E-2 

2 

2.104E3 

2.583E3 

2.279E1 

3 

3.506E3 

1.390E4 

2.964E2 

Table 3. 

Cantilever Frequency Comparisons (n = 2) 

Mode 

Theoretical 

Computed 

Error 

Number 

Frequency (Hz) 

Frequency (Hz) 

% 

1 

7.012E2 

7.011E2 

-2.169E-2 

2 

2.104E3 

2.239E3 

6.440E0 

3 

3.506E3 

7.043E3 

1.009E2 

Table 4. Cantilever Frequency Comparisons (n = 3) 

Mode 

Theoretical 

Computed 

Error 

Number 

Frequency (Hz) 

Frequency (Hz) 

% 

1 

7.012E2 

7.012E2 

-2.170E-2 

2 

2.104E3 

2.167E3 

3.003E0 

3 

3.506E3 

3.879E3 

1.063E1 
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Table 6. Cantilever Frequency Comparisons (n = 4) 


Mode 

Number 

Theoretical 
Frequency (Hr) 

Computed 
Frequency (Hz) 

Error 

% 

1 

7.012E2 

7.010E2 

-3.235E-2 

2 

2.104E3 

2.120E3 

7.937E-1 

3 

3.506E3 

3.680E3 

4.950E0 


The accuracy of the computed frequencies is improved by taking into account the 
higher order correction terms. However, when n ^ 6, the reduced mass matrix is no longer 
positive definite and the eigenvalue solution process breaks down. This limitation of adding 
a finite number of correction terms can be explained by the fact that the third frequency 
of the overall structure exceeds the lowest frequency of the o mi t get (O-set) system thus 
violating the convergence criterion for Guyan reduction. 

Another cantilever example is constructed by assuming that the three elements near 
the free end are made up of a material with E = 30 x 10 4 psi instead of steel. By retaining 
the same degrees of freedom as in the previous example of all steel construction, it becomes 
possible to add an almost limitless number of correction terms. This is because there is 
no overlap between the frequency spectrum of the first three modes of the full system and 
that of the O-set system. 

6.2 Plama Motor Generator ( PMG): This example comes from the modal testing 
and finite element analysis of the PMG Far End Package (Figure 1). The PMG experiment 
is a payload on a Delta II 7925 launch vehicle. The mission is scheduled to take place in 
July 1993. 



Figure 1. PMG Far End Package 
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The analysis set degrees of freedom correspond to the accelerometer locations used in 
the modal survey test. The improved Guyan reduction is performed with different n on 
the PMG Far End Package model. The computed frequencies are compared with those of 
the full model which are taken as the reference solution and the results are listed in Tables 
6 through 8. Several frequencies that were not found by the standard Guyan reduction 
start to reappear by adding the correction terms. 


Tfeble 0. PMG Frequency Comparisons (n = 0) 


Mode 

Reference 

Computed 

Error 

Number 

Frequency (Hz) 

Frequency (Hz) 

% 

1 

56.32 

56.39 

H 

2 

84.46 

84.66 

ro 

3 



wM 

4 

118.19 

119.58 

1.17 

5 

159.46 

160.48 

0.63 

6 


— 


7 

185.19 

— 


8 

215.16 

220.65 


9 

217.65 

224.09 


10 

228.36 

236.73 

3.56 

11 

234.52 

— 

— 

12 

243.43 

256.55 

5.38 

13 

264.53 

— 

— 

14 

299.03 

— 

— 

15 

305.16 

330.53 

8.31 

Table 7. PMG Frequency Comparisons (n 

= 1) 

Mode 

Reference 

Computed 

Error 

Number 

Frequency (Hz) 

Frequency (Hz) 

% 

1 

56.32 

56.32 

0.00 

2 

84.46 

84.46 

0.45E-4 

3 

100.66 

100.66 

0.59E-3 

4 

118.19 

118.20 

0.46E-2 

5 

159.46 

159.46 

0.16E-2 

6 

170.06 

171.65 

0.92 
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Table 8. PMG Frequency Comparisons (n = 2) 


Mode 

Number 

Reference 
Frequency (Hz) 

Computed 
Frequency (Hz) 

Error 

% 

1 

56.32 

56.82 

0.89 

2 

84.46 

85.52 

1.24 

3 

100.66 

101.12 

0.45 

4 

118.19 

118.31 

0.09 

5 

159.46 

160.32 

0.53 

6 

170.06 

170.18 

0.064 

7 

185.19 

185.29 

0.054 

8 

215.16 

215.19 

0.016 

9 

217.65 

217.71 

0.026 

10 

228.36 

228.61 

0.10 

11 

234.52 

234.60 

0.035 

12 

243.43 

242.79 . 

-0.26 

13 

264.53 

264.65 

0.043 

14 

299.03 

299.59 

0.17 

15 

305.16 

305.60 

0.14 


6. Conclusion: 

A noniterative procedure to enhance the standard Guyan reduction with a series of terms 
has been presented. In practice, it may be possible to add only a finite number of the 
correction terms as demonstrated by the NASTRAN examples. 
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